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摘 要 建立 了 Ti6AI4V 合 金 铸 件 / 铸 型 界面 换 热 系数 ( 门 的 一 维 反 算 模型 , 从 数学 及 数值 模拟 的 角度 研究 了 型 壳 热 物性 参数 
和 热电 偶 定位 等 参数 对 hh 计 算 的 影响 , 分 析 了 不 同 参数 影响 的 不 同 特点 , 据 此 对 型 壳 热 物 性 参数 和 热电 偶 定位 位 置 等 进行 
了 修正 , 提高 了 hn 反 算 精度 . 修正 计算 参数 后 的 反 算 结果 表明 , Ti6Al4V 合金 熔 模 铸造 条 件 下 ,天 的 变化 可 分 为 4 个 阶段 : (1) 
铸件 为 液态 ,维持 约 440 W/m KJ); (2) 铸件 表面 生成 完整 凝固 层 , 此 阶段 4 下 降 近 60%; (3) 凝固 层 不 断 增 厚 至 铸件 凝固 ， 
此 阶段 A 下 降 接 近 峰 值 的 20%; (4) 铸件 凝固 后 , 产 随 温度 缓慢 下 降 . 在 三 维 模型 中 对 反 算 得 到 的 /进行 了 验证 , 得 到 的 模拟 
温度 与 实测 温度 基本 吻合 , 表明 反 算 得 到 的 /较为 准确 , 可 以 应 用 于 Ti6Al4V 合金 熔 模 铸造 过 程 的 数值 模拟 中 ， 
关键 词 钛 合金 , 熔 模 铸造 , 界面 换 热 系数 

中 图 法 分 类 号 TG249.5 文献 标识 码 A 文章 编号 0412-1961(2015)08-0976-09 


RESEARCH ON THE INTERFACIAL HEAT TRANSFER 
COEFFECIENT BETWEEN CASTING AND CERAMIC 
SHELL IN INVESTMENT CASTING PROCESS 

OF Ti6Al4V ALLOY 


SHAO Heng », LI Yan », NAN Hai™, XU Oingyan " 
1) Key Laboratory for Advanced Materials Processing Technology, Ministry of Education, School of Materials Sci- 
ence and Engineering, Tsinghua University, Beijing 100084 

2) Beijing Institute of Aeronautical Materials, Beijing 100095 

Correspondent: XU Qingyan, professor, Tel: (010)62795482, Email: scjxqgy@mail.tsinghua.edu.cn 

Supported by China- EU (European Union) Science & Technology Cooperation in Aviation, 7th 
Framework Programme of EU, National Basic Research Program of China 
(No.2011CB706801), National Natural Science Foundation of China 
(Nos.51171089 and $51374137), National Science and Technology Major Project 
(No.2012ZX04012011) and High Technology Research and Development Program 
of China (No.2007AAO04Z141) 

Manuscript received 2015-01-12, in revised form 2015-05-02 


ABSTRACT Investment casting process is an important way to get complex parts of titanium alloy. However 
there are few research on the interfacial heat transfer coefficient (h) between casting and shell thus the temperature 


simulation of investment casting process of titanium alloy is often inaccurate. In order to get a relatively accurate h, 
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钛 合金 具有 比 强 度 高 、 耐 热 性 好 、 耐 刨 性 好 、 生 ” 距 的 影响 , 并 反 算 得 到 各 工艺 参数 (型 壳 预 热 温 度 和 
物 相 容 性 好 等 优点 , 因而 被 广泛 应 用 于 航空 航天 、 ”型 这 厚度 ) 下 铸件 凝固 过 程 中 界面 换 热 系数 随 铸件 


其 邵 瑜 等 : 熔 模 铸造 条 件 下 Ti6Al4V 合金 铸件 与 陶瓷 型 沉 间 界面 换 热 系 数 研究 


a one-dimensional mathematical model for the reverse calculation of h between casting and shell in investment cast- 
ing process of Ti6Al4V alloy was built and the analytic relationship between temperature and time/heat flux was es- 
tablished. Considering the calculated h is significantly affected by the error of parameters such as the specific heat 
capacity and thermal diffusivity of shell and position of thermocouples, research on the error of these parameters is 
essential. The relationship between the error of these parameters and the temperatures in the casting and shell was 
studied and it was found that the effect of different kind of error on the temperature field was obviously different. 
An experiment based on the one-dimensional mathematical model was done and temperatures of different positions 
were measured. Based on the effect of different kind of error and the difference between the calculate temperature 
field and the measured temperatures, the proportion of effect of each kind of error was assessed. These errors were 
revised on the basis of the assessment, thus a relatively accurate h between the casting and shell was obtained. The 
relationships between h and thickness of the solidified layer on the casting/temperature at the surface of casting can 
be divided into 4 stages: (1) Metal was liquid and h kept about 440 W/(m’* K); (2) Solid layer appeared on the sur- 
face, and h declined nearly 60%; (3) Solid layer grew up before metal became completely solid and h declined near- 
ly 20% of its maximum; (4) After metal solidified, h declined slowly as temperature on the surface of casting 
dropped. These relationships were applied in a three-dimensional model for numerical simulation of the tempera- 
ture field. Temperatures of different positions in casting and shell were calculated and calculated temperatures 
agreed with measured ones well. Thus the accuracy of h was identified and it can help solve problems in the pro- 
duction in investment casting process of Ti6Al4V alloy. 


KEY WORDS titanium alloy, investment casting, interfacial heat transfer coefficient 
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医疗 及 化 学 工业 . 但 钰 合金 零件 往往 具有 十 分 复杂 ” 温度 的 变化 . Sun 等 "研究 了 挤 压铸 造 条 件 下 


到 


但 


造 工 艺 来 成 形 . 在 众 模 铸造 工艺 开发 过 程 中 , 数值 模 


时 间 和 成 本 而 受到 青睐 , 得 到 越 来 越 多 的 应 用 


四 


的 形状 , 难于 进行 机 械 加 工 , 因而 需要 采用 熔 模 铸 AM60 镁 合金 与 金属 铸 型 之 间 的 界面 换 热 行为 , 比 
较 多 项 式 拟 合 方法 与 反 算法 计算 得 到 的 界面 换 热 
因为 能 够 预测 缺陷 的 形成 , 减少 工艺 开发 所 需 。 系数 , 结果 表明 反 算 法 的 精度 较 高 . 徐 宏 等 ”研究 了 
.和 砂 型 铸造 条 件 下 ZL104 与 铸 型 间 的 界面 换 热 系数 ， 
数值 模拟 的 精确 度 容易 受到 相关 模拟 参数 精确 。 给 出 了 界面 换 热 系数 随时 间 的 变化 . Wang 等 "测量 


度 的 直接 影响 , 如 铸件 和 铸 型 的 热 物 性 参数 .铸件 了 连续 铸造 条 件 下 到 弯 月 面 不 同 距离 处 的 钢材 表 


铸 型 间 的 界面 换 热 系数 等 . 通常 而 言 , 前 者 可 以 通 


挤 


口 
口 


一 


何 


查 | 


六 文献 或 实验 测 和 


钛 合金 熔点 高 , 化 学 活性 高 , 多 采用 熔 模 铸造 


面 温 度 , 反 算 得 到 连续 铸造 过 程 中 的 界面 换 热 系 
EE 得 到 较为 准确 的 数值 , 后 者 。 数 , 并 以 此 预测 不 同位 置 钢材 表面 凝 壳 厚度 , 预测 
都 是 通过 估计 给 出 , 且 经 常 作 为 常数 使 用 "可靠 。 结果 与 实测 结果 吻合 良好 . 
性 和 精确 性 都 有 待 商 枞 . 因此 , 对 熔 模 铸造 条 件 下 
钛 合金 铸件 与 铸 型 界面 换 热 系 数 的 专门 研究 就 最 ”工艺 , 铸 型 采用 陶瓷 材料 , 铸造 过 程 在 真空 炉 中 完 
得 尤为 必要 . 成 ,大 大 增加 了 实验 设计 与 数据 处 理 的 难度 , 目 


近年 来 , 国内 外 学 者 对 铸造 过 程 铸件 / 铸 型 界面 ”前 尚 鲜 有 对 其 铸造 过 程 铸件 / 铸 型 界面 换 热 系数 的 


换 热 系数 的 研究 , 主要 涉及 压铸 后 、 熔 模 铸 造 "3、 研究 , 无 法 为 生产 过 程 提 供 更 多 的 帮助 . 因此 , 本 工 
压铸 造 " 和 砂 型 铸造 "和 连续 铸造 ?2 等 , 包含 作 拟 采用 热传导 反 算 法 , 测算 Ti6Al4V 合金 熔 模 铸 


合金 ,位 合金 .高温 合金 和 钢 等 , 取得 了 较 好 的 效 ” 造 铸件 /型 壳 界 面 换 热 系数 , 以 期 为 Ti6Al4V 合金 熔 


果 . 郭志 脑 等 "针对 铝 合 金 压铸 条 件 下 的 界面 换 热 。 ” 模 铸 造 过 程 数值 模拟 提供 准确 的 计算 参数 . 
为 展开 了 大 量 的 研究 , 建立 了 完善 的 数值 模型 ， 


1 物理 /数学 模型 


分 析 了 铸件 厚度 、 压 铸 压 力 、 铸 型 温度 和 测量 位 置 

等 对 界面 换 热 系数 的 影响 . Konrad 等 "针对 ”1.1 铸件 / 铸 型 界面 换 热 物理 模型 
IN738LC 高温 合金 熔 模 铸造 条 件 下 的 界面 换 热 行 本 工作 采用 平板 形 铸件 , 测量 平板 中 央 部 位 型 
为 , 测量 了 陶瓷 型 壳 的 比热容 和 热 扩散 率 , 研究 了 ” 吝 的 温度 . 考虑 到 铸件 高 度 与 宽度 远大 于 厚度 , 因此 
型 壳 预 热 温 度 和 厚度 对 铸件 晶 粒 尺寸 和 二 次 臂 间 ”将 平板 中 央 部 位 型 过 的 温度 场 近似 视 为 随 距 铸 件 / 型 
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过 界面 的 距离 而 变化 的 一 维 温度 场 , 如 图 1 所 示 . 
在 型 壳 中 , 其 温度 7 分 布 为 : 


T=T(, 7, gg (1) 
q1= 47) (2 
gq = 92(7) G3) 


图 1 及 式 (1)~(3) 中 ,7 为 型 过 中 不 同位 置 、 不 同时 刻 
的 温度 ,x 为 到 铸件 /型 这 界面 的 距离 , z 为 时 间 , qi 
为 铸件 /型 这 界面 的 热流 密度 , 9 为 型 壳 外 表面 热 
1.2 铸 型 温度 分 布 数 学 模型 

在 简化 后 的 一 维 传 热 模 型 中 , 铸 型 内 的 温度 分 
布 满足 以 下 条 件 : 


oo7_ 39,,07 
| 全 
T(x, 0， di, 0) 二 1 (5) 
aT 
AZ o= -qi (O 


后 | ed 
”一 上 o oo 
， Casting Mold Vaccum 
0o 0 oo 
OO | 
ms 省 
0o oO | 
oo oo © 一 
oo © © 
一 
0 d ” 
(b) 
T 
Casting Vaccum 


图 1 一 维 传 热 模型 及 其 温度 场 示 意图 


Fig.1 Schematics of one- dimensional heat transfer model 
(a) and temperature field in the model (b) (9 一 inter- 
facial heat flux between casting and mold; 9g: 一 heat 
flux across the surface of mold; x—distance from 
the interface between casting and mold; 4d—thick- 


ness of mold; 7 一 temperature) 
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AS 全 -dg2(7T) (7) 


式 (4)~(7) 中 ，p, ce 和 入 分 别 为 型 这 材料 的 密度 、 比 热 
容 和 导热 系数 ; 7 为 型 壳 初 始 温 度 . 

对 陶瓷 型 沉 热 物性 参数 的 实际 测量 结果 表明 ， 
型 壳 的 c 与 入 随 温度 的 变化 均 较 小 . 忽略 温度 对 型 
壳 热 物性 参数 的 影响 , 将 p,c 和 A 作为 常数 处 理 , 式 
(4) 可 简化 为 : 


pe -5 去 (8) 


暂 令 gq,=0, 仅 考 虑 0 对 7 的 影响 , 则 有 : 


T=T(x, 7, q) (9) 

T(x, 0, q)= 7 (10) 
oT 

A l=0 (11) 


式 (9) 中 的 了 满足 式 (90), (8), (10) 和 (11) 的 各 项 条 件 . 

假设 对 某 一 特定 的 热流 g', 当 gq =gi 时 ， 
T=7" 为 满足 式 (6), (8), (10) 和 (11) 的 解 . 将 9” 任意 
分 拆 为 2 股 热流 gi! 和 g?, 即 在 任意 时 刻 均 满足 
gi =gu+g7 ,假设 T=7T” 和 7=7 分 别 为 gq,=gl 和 
qi=gi 时 满足 式 (6), (8), (10) 和 (11) 的 解 , 令 
7 =7*+7 -7 , 则 7 满足: 


X90= -9 (7) -gr (7) = (12) 


Tx, 0, q)=T,+7, -7,=7, (13) 
5 
As =0+0=0 (14) 


由 于 式 (8) 中 p,c 和 和 作为 与 温度 .时 间 无 关 的 
常数 处 理 , 故 7， 也 满足 式 (8). 式 (12) 即 为 w = 时 
的 式 (6), 而 式 (13) 和 (14) 与 式 (10) 和 (11) 相 同 , 即 7 
也 为 mw = 时 满足 式 (6), (8), (10) 和 (11) 的 解 , 由 解 
的 唯一 性 知 T=T .由 以 上 可 知 , 当 q =gi+gr 
时 , gi 对 温度 场 造成 的 影响 (7" -7) 也 相当 于 gi 和 
gr 对 温度 场 造成 的 影响 (7-7) 和 (7, -7 的 侍 
加 , 同 理 , 将 wo 看 做 MM 股 (M 为 任意 正 整 数 ) 热 流 的 


M 
合 加 》 gw” 时, 即 : 
m=1 


M 
pA 而 (15) 
才 % 窗 量 


若 7=7, 是 m=g" 时 满足 式 (6), (8), (10) 和 
(11) 的 解 , 则 有 : 


M 
7 -7 = 7, -70 (16) 
m=1 
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M 
亦 即 若 将 gi 看 做 MM 股 热流 的 难 加 了 gq”, 则 这 
m=1 


些 热流 对 温度 场 的 影响 相互 独立 , 总 的 效果 为 各 热 
流 的 效果 的 简单 合 加 . 同样 的 , 对 于 g, 而 言 , 若 将 q， 
看 做 若干 股 热流 的 县 加 , 则 这 些 热流 对 温度 场 的 影 
响 也 相当 于 各 热流 的 影响 的 简单 合 加 , 而 且 0 和 0， 
对 于 型 这 温度 场 的 影响 相当 于 和 9 各 自 对 温度 
场 的 影响 的 简单 县 加 . 
1.3 表面 换 热 系数 求解 方法 

界面 换 热 系数 的 计算 采用 Beck 等 所 提出 的 非 
线性 估算 法 , 其 基本 原理 是 : 在 某 一 时 刻 r 之 前 界 
面 换 热 系数 已 经 计算 完毕 (已 知 ); 假设 从 这 一 时 刻 
起 NW 个 时 间 步 长 Ndt 内 界面 换 热 系数 保持 不 变 , 
其 值 为 ,基于 之 前 的 界面 换 热 系数 和 温度 场 , 以 
及 此 时 的 界面 换 热 系数 及 计算 未 来 N 个 时 间 步 长 
内 温度 场 的 变化 , 即 : 

i (1 


T+Ndr 


T=T(x, 7’ +dr7, gog] (18) 


式 (19) 中 = 六 CT -70) =ce(7 -75')( 其 
中 , 7,, 为 计算 得 到 的 铸件 表面 温度 ; o 和 s 为 有 
效 辐射 率 和 黑体 辐射 常数 ; 7 为 背景 温度 , 即 真 空 
炉 内 壁 温度 )， 

定义 f(h) 为 : 


MK 
/(n) = > Toor + nd7, qi, 9;)— 7 (19) 
nhl 


式 (19) 中 ,7” -ar 为 了 +tndr 时 刻 % 点 测 得 的 温 
度 . jj) 代表 二 时 刻 之 后 界面 换 热 系数 保持 为 六 
的 情况 下 计算 得 到 的 温度 场 与 测量 所 知 的 温度 场 
之 间 的 差异 . /Fe) 越 小 , 则 表明 对 二 时 刻 之 后 的 界 
面 换 热 系 数 估计 值 h 越 接近 真实 情况 . 因此 计算 过 
程 以 极 小 化 AU) 为 目的 , 当 f(4) 达 到 极 小 值 时 的 
六 即 为 时 刻 所 求 的 界面 换 热 系 数 ， 
2 实验 设计 
本 工作 采用 了 左右 对 称 的 2 块 板 状 铸件 , 以 利 
于 获取 较 多 实验 数据 . 实验 采用 的 模型 及 热电 偶 位 
置 安排 如 图 2 所 示 . 铸件 板 状 部 分 尺寸 为 170 mmx 
170 mmx25 mm. 不 同 热电 偶 分 别 埋 入 铸件 或 型 这 
中 . 埋 入 铸件 的 热电 偶 , 编号 从 上 到 下 依次 为 TC 
casting 1, 2, 3 和 4, 分 别 位 于 冒 口 ` 板 状 铸件 中 央 、 
内 浇 道 和 直 浇 道中 , 其 中 TC casting 1, 2 和 4 号 热 
电 偶 采 用 钨 铁 热 电 偶 , 用 来 测量 铸件 温度 变化 , TC 


casting 3 号 热电 偶 为 K 型 热电 偶 , 用 来 测量 金属 液 
到 达 时 间 . 埋 入 型 壳 中 的 热电 偶 处 于 同一 水 平 高 
度 , 横向 间距 10 mm, 位 于 板 状 铸件 中 央 附 近 , 热电 
偶 顶 端 到 铸件 /型 这 界面 距离 依次 为 3, 6 和 9 mm， 
按照 离 界 面 由 近 到 远 其 编号 依次 为 TC mold 1, 2 和 
3, 均 为 K 型 热电 侦 . 铸件 本 身 为 对 称 造型 , 在 对 侧 
于 同样 位 置 放置 同样 的 热电 偶 , 可 在 一 次 浇注 中 获 
得 2 组 数据 用 于 反 算 . 

实验 采用 的 型 过 材料 热 物性 通过 实验 测 得 . 型 
过 材料 比热容 随 温度 的 变化 如 图 3 所 示 , 其 在 
900 左右 的 突然 上 升 与 型 过 中 高 岭 土 成 分 的 莫 来 
石化 有 关 . 型 过 材料 导热 系数 受 温度 影响 很 小 , 基 
本 保持 在 1.4 W/(m.K) 左 右 , 密度 则 为 2170 kg/m. 
实验 采用 的 Ti6Al4V 合金 名 义 成 分 (质量 分 数 ) 为 : 
Al 6%, V 4%, 余 量 为 Ti, 其 比热容 及 导热 系数 随 温 
度 的 变化 见 图 4. Ti6Al4V 比热容 和 导热 系数 在 
1620~1650 CC 的 突变 对 应 于 熔化 /凝固 过 程 , 而 980~ 
1000 C 的 突变 对 应 于 固态 相 变 . 


2 铸 型 型 腔 形状 及 热电 偶 布 置 示意 图 
Fig.2 Schematic of the cavity and the positions of the ther- 


mocouples (TC) 
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3 型 壳 材 料 比热容 (c) 随 温度 的 变化 
Fig.3 Variation of specific heat capacity (c) of shell with 
temperatures (7) 


由 人 
«Iz 


ChinaXiv 合 作 期 刊 


980 金 属 学 报 第 51 卷 
900 
(a) ee 1600| 六 
800 上 | ' lc casting 4 
a | 一 TC casting 1 
a OF 人 i 1200 上 | TC casting2 
立 Am i | : 
“bp pa ?2 | 
600 上 i SS | 
兰 于 ~ 800 | 
> 500 上 | 
400[ | 
400[ | 
300 ! 1 ! 0 1 | 1 | 
0 500 1000 1500 2000 0 1000 2000 3000 4000 
ya / °C 过 入 
图 5 实验 测 得 铸 型 及 铸件 各 处 的 温度 
4 Fig.S Measured temperatures at different positions in 
四 mold and casting (+—time) 
30 上 | | 
网 有 
国 一 
一 300 上 300 
三 20 4 
三 , 兰 200 
~ ee 加 § 
-一重 - 
10| mm 站 “200 三 100 
CE 三 并 
三 0 
0 | 1 | 0 50 10 150 200 
0 500 1000 1500 2000 100r r/s 
T/°C 
图 4 Ti6AI4V 合 金 的 c 和 导热 系数 4 随 温度 的 变化 Wo 100 2000 3000 4000 
Fig.4 Variations of c (a) and thermal conductivity 4 (b) of r/s 


™ 
> 
CN 
Co 
LO 
©O 
© 
CO 
3 
可 
© 
CN 


Ti6Al4V with temperatures 


6 反 求 得 到 的 界面 换 热 系 数 (有 ) 及 局 部 放大 图 


| 
[4 


chinaXiv 


Fig.6 Calculated interfacial heat transfer coefficient (1) be- 
tween 0~4000 s and between 0~200 s (inset) 


3 结果 与 讨论 
3.1 实验 测 温 及 反 求 结果 
图 5 是 实验 测 得 的 一 侧 热电 偶 的 温度 . 根据 多 ”的 各 处 温度 变化 , 结合 铸件 、 铸 型 的 热 物 性 参数 和 
次 测量 的 数据 , 热电 偶 不 带 保 护 套 时 浇注 后 数秒 内 初始 状态 , 各 热电 偶 的 定位 位 置 , 可 以 通过 反 算 精 
温度 升 至 1700 'C 左 右 , 而 带 保护 套 时 升温 较 慢 ,最 。 确 地 得 到 h 的 变化 . 然而 实际 情况 下 铸 型 的 热 物性 
高 温度 在 约 1570 'C 处 停留 较 长 时 间 , 而 Ti6AIM4V 合 ” 参数 及 热电 偶 的 位 置 并 不 总 是 非常 准确 , 而 OIMa- 
金 两 相 区 为 1620~1650 'C, 取 其 中 值 1635 'C, 估计 honey 等 指出 , 这 些 误差 对 反 算 结 果 会 造成 很 大 的 
出 2 号 热电 偶 处 测量 得 到 的 温度 大 约 比 实际 温度 低 。 ”影响 , 因此 有 必要 评估 这 些 误差 造成 了 多 大 的 影响 
4% 左 右 . 而 型 壳 中 的 热电 偶 工作 端 均 随同 型 壳 一 起 。 并 尽量 消除 它们 的 影响 . 
烧结 , 与 型 这 接触 较 好 , 因此 认为 测 得 的 温度 即 为 图 7~9 分 别 是 当 计 算 采 用 的 铸 型 材料 c、 热 扩 
相应 位 置 型 壳 温 度 . 图 6 是 根据 TC casting2 和 TC ” 散 率 a 和 热电 偶 实 际 位 置 偏离 预定 位 置 时 计算 得 
mold 1~3 号 热电 侦 的 温度 反 求 得 到 的 h-t 关 系 图 ， ”到 的 与 预 设 h 之 间 的 关系 以 及 在 条 件 A, C,E 下 
其 中 插图 是 对 h-t 关 系 在 0~200 s 内 变化 情况 的 局 。” 反 算 过 程 中 得 到 的 铸件 、 铸 型 温度 与 此 前 无 误差 情 


部 放大 . 况 下 模拟 得 到 的 温度 的 对 比 . 图 中 条 件 A~F 依次 
3.2 界面 换 热 系数 反 求 精度 评估 为 : 计算 采用 的 c 偏 小 和 偏 大 20%, 计算 采用 的 w 
3.2.1 铸 型 材料 比热容 、 热 扩散 率 及 热电 偶 定 位 偏 小 和 偏 大 20%, 以 及 实际 的 热电 偶 位 置 为 2.5， 


误差 对 计算 结果 影响 评估 在 一 维 传 热 模 型 中 ，, 
知 铸 件 / 铸 型 之 间 产 随 铸件 状态 的 变化 关系 已 知 , 铸 
件 、 铸 型 的 热 物性 参数 准确 且 初 始 状态 已 知 , 各 热 
电 侦 定位 准确 , 测 得 的 温度 无 误差 , 通过 温度 场 模 
拟 可 以 得 到 各 热电 偶 处 的 温度 变化 . 根据 模拟 得 到 


5.5, 8.5 mm 和 3.5, 6.5, 9.5 mm. 

对 比 图 7~9 表 明 , 铸 型 材料 c a 和 热电 侦 实 际 
位 置 偏离 预定 位 置 带 来 的 的 误差 和 反 算得 到 的 温 
度 变化 误差 有 明显 不 同 的 特点 : c 的 误差 对 铸件 温 
度 和 铸 型 中 各 点 的 温差 有 明显 的 影响 , 而 且 会 对 天 
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9 热电 侦 埋 入 型 沉 位 置 发 生 偏离 时 反 算得 到 的 有 与 模拟 温度 
Fig.9 Calculated h based on conditions E and F (when real distances between the interface and thermocouples are 2.5, 5.5, 
8.5 mm, and 3.5, 6.5, 9.5 mm) (a) and temperatures based on condition E (b) 
的 计算 造成 很 大 的 误差; a 的 误差 对 铸件 温度 影响 差 基 本 没有 影响 . 
不 大 , 而 对 铸 型 中 各 点 温差 有 明显 的 影响 , 对 h 的 计 3.2.2 误差 评估 与 校正 基于 对 型 过 材料 
算 精 度 影 响 不 大 ; 热电 偶 定 位 误差 对 铸件 温度 和 h ” 比热容 和 导热 系数 的 测量 结果 , 可 认为 c=1290 J/ 


的 计算 精度 的 影响 几 


乎 完全 体现 在 浇注 开始 后 约 


50 s 内 , 对 后 续 反 算 的 影响 也 完全 继承 自 最 初 30 s 


内 的 影响 , 同时 热电 偶 定 位 误差 对 铸 型 内 各 点 的 温 


(kg:K), 4=1.4 W/m.K), ao=5.0x107m2s. 实验 中 型 


壳 平 均 厚 度 约 为 18 mm, 热 


EE 侦 与 界面 间 的 等 


效 最 


小 距离 为 3.5 mm, 浇注 温度 约 为 1720 'C( 考 虑 到 
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Ti6Al4V 合金 液 进 入 型 腔 之 前 在 直 浇 道 和 横 浇 道 中 
与 铸 型 接触 , 同时 本 实验 采用 重力 浇注 , 流速 较 慢 ， 
认为 Ti6Al4V 合金 进入 型 腔 后 的 平均 温度 约 为 
1680 °C). 

图 10 为 0~750s 内 及 反 求 过 程 中 计算 得 到 的 各 
测 温 点 温度 与 实测 温度 的 对 比 . 模拟 得 到 的 型 过 最 
高 温度 约 为 840 'C, 可 以 认为 型 过 中 高 岭 土 成 分 的 
莫 来 石化 影响 很 小 . 从 图 10 可 以 看 出 , 按照 原 有 的 
型 壳 热 物性 参数 计算 得 到 的 型 沉 中 不 同位 置 热 电 
偶 的 温差 明显 小 于 实测 不 同 热电 侦 的 温差 , 铸件 温 
度 下 降 速 率 明 显 低 于 实际 值 , 从 之 前 的 评估 可 以 认 
为 采用 的 型 过 w 高 于 实际 值 , 而 c 低 于 实际 情况 . 此 
外 , 铸件 温度 开始 明显 下 降 的 模拟 时 间 明 显 低 于 实 
际 时 间 , 表明 热电 偶 到 界面 的 最 小 等 效 距 离 (与 距 界 
面 最 近 的 热电 侦 温度 相等 的 型 这 内 部 温度 点 到 界 
面 的 距离 ) 低 于 实际 值 . 

针对 上 述 问 题 , 将 c 适 当 上 调 , 将 a 适当 下 调 ， 
同时 微调 热电 偶 到 界面 的 最 小 等 效 距离 以 保证 反 
算得 到 的 铸件 温度 与 实测 值 相 符 以 及 铸 型 中 不 同 
测 温 点 温差 与 实测 温差 相符 . 经 多 次 调整 后 , 当 c= 
1410 J/(kg-K), o=:4.09x107mzs, 热电 偶 到 界面 的 最 
小 等 效 距离 为 3.6 mm 时 ， 反 算得 到 的 铸件 钛 型 温度 
与 实测 值 吻合 较 好 , 图 11 为 修正 后 反 算得 到 的 hh 及 
铸件 铸 型 的 温度 . 

从 图 11 可 以 看 出 , 修正 铸 型 材料 热 物 性 参数 
及 热电 偶 定 位 位 置 后 反 算 过 程 中 得 到 的 铸件 及 铸 
型 温度 与 实测 温度 吻合 良好 . 从 反 算 得 到 的 有 随 
时 间 变 化 的 关系 来 看 , 浇注 开始 后 约 20 s 内 瑚 变 
化 剧烈 , 从 0 上 升 至 约 440 W/m .BR), 很 快 又 下 降 
至 约 184 W/m.K), 对 应 于 金属 液 开 始 与 铸 型 接触 
并 于 表面 生成 一 层 凝 壳 . 随后 随 着 凝固 的 进行 ,组 
慢 下 降 至 约 110 W/m .KR), 这 一 阶段 随 着 凝 壳 的 增 
厚 和 铸件 的 收缩 , 凝 壳 与 铸 型 间 的 接触 状况 不 断 恶 
化 ,天 的 变化 量 相 当 于 峰值 的 20%. 铸件 完全 凝固 
后 ,h 随 着 铸件 的 冷却 而 降低 , 这 一 阶段 由 于 铸件 不 
有 明显 的 收缩 ,hn 的 下 降 速 度 较 前 一 阶段 还 要 低 
得 多 . 

3.3 界面 换 热 系数 分 析 

3.3.1 界面 换 热 系数 与 铸件 状态 的 关系 如 

上 节 所 述 , 浇注 后 , 铸件 状态 的 变化 经 过 以 下 几 个 
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10 反 求 得 到 的 各 测 温 点 温度 与 实测 温度 对 比 


Fig.10 Comparison between measured and calculated tem- 


peratures 
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11 修正 铸 型 材料 c  w 和 热电 偶 定位 位 置 之 后 反 求 得 

到 的 及 和 各 测 温 点 温度 与 实测 温度 对 比 

Fig.11 Calculated h between 0~4000 s and between 0~200 s 


(inset) (a), comparison between measured and calcu- 


lated temperatures with specific heat capacity, ther- 
mal diffusivity of shell and positions of thermocou- 


ples adjusted (b) 


相 线 的 部 分 为 凝固 层 , 则 可 以 得 到 铸件 表面 凝固 层 
量度 5 随时 间 的 变化 . 

基于 h 和 6 随时 间 的 变化 , 可 以 得 到 这 2 者 之 间 
的 关系 , 如 图 12 所 示 . 图 12 表明 , 铸件 表面 形成 凝 


阶段 : 完全 液态 , 表面 生成 凝固 壳 层 , 凝固 层 厚 度 不 
断 增 大 , 铸件 完全 凝固 . 反 算 不 仅 可 以 得 到 有 随时 间 
的 变化 , 也 可 以 得 到 铸件 内 部 的 温度 、 凝 固 分 数 分 
布 . 在 铸件 完全 凝固 前 , 以 温度 低 于 Ti6Al4V 合金 液 


回 层 前 ,h 仍 处 于 峰值 附近 , 而 凝固 层 一 旦 形成 , 则 有 
急剧 下 降 . 当 凝 固 层 ee oe 有 降 至 约 
184 W/m .KR), 随后 随 着 凝固 层 厚度 的 增加 而 缓慢 

下 降 . 当 凝 固 层 厚度 达到 11 ne 1 与 完全 凝固 时 
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图 13 铸件 基本 凝固 后 4 随 铸件 温度 变化 的 关系 
Fig.13 Relationship between h and T of the casting after 
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形成 , 铸件 与 型 充 


和 2 段 中 , 金 
壳 , 糊 状 区 向 铸件 内 部 


属 液 与 铸 型 的 表面 开始 形成 凝 
进 . 由 于 完整 的 凝 壳 逐 渐 
的 接触 状态 开始 发 生 改 变 , 因此 及 


[ES 


14 根据 反 求 得 到 的 疡 计算 出 的 型 过 和 铸件 温度 
Fig.14 Simulated temperatures of casting and shell based 
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急剧 减 小 , 降幅 近 60%. 第 3 段 中 , 铸件 表层 已 经 形 
成 完整 的 凝 壳 , 其 与 型 壳 的 接触 状态 固定 下 来 , 不 
再 发 生 明 显 改变 , 因此 降低 的 速率 大 大 减 慢 , 这 一 ”测量 到 的 温度 基本 相符 , 表明 反 求 得 到 的 刀具 有 较 
阶段 4 的 降幅 约 为 其 峰值 的 20%. 此 后 , 铸件 基本 凝 。 好 的 准确 度 , 也 证 明 选 取 的 计算 参数 具有 较 高 的 准 
固 , 温度 开始 快速 下 降 , 此 后 有 主要 取决 于 界面 两 侧 确 度 . 同时 , 计算 得 到 的 温度 与 测量 到 的 温度 间 有 
的 温度 和 辐射 率 等 . 部 分 位 置 有 一 定 差异 , 表明 本 反 求 模型 仍然 存在 一 

铸件 凝固 后 , 基于 hh 和 铸件 表面 温度 7 随时 间 ” 定 的 误差 . 造成 这 些 误差 的 原因 主要 有 : (1) 实验 中 
的 变化 , 可 以 得 到 它们 之 间 的 关系 , 如 图 13 所 示 . 图 型 系统 并 非 标准 的 一 维 模型 , 从 而 使 反 算 存 
13 表明, 铸件 基本 凝固 之 后 , 其 与 型 沉 间 的 h 随 铸件 在 系统 生 误 差 ; (2) 型 碗 与 铸 型 本 身 形 状 的 不 精确 , 
温度 下 降 , 其 间 昌 有 波动 , 但 影响 不 大 . 铸件 凝固 后 。 以 及 布置 的 热电 侦 对 温度 场 的 影响 , 造成 温度 偏离 
7 与 铸件 温度 间 的 关系 可 简化 为 : 实际 情况 ; (3) 型 壳 、 铸 件 热 物性 参数 及 热电 偶 位 置 

h=4.2584x 1037s -6.5324x 103572 + 等 实验 参数 的 偏差 , 从 而 导致 计算 与 实际 不 一 致 . 


4 结论 


出 , 根据 反 求 得 到 的 计算 出 的 型 这 和 铸件 温度 与 


5.0115 x 1027+21.478 (21) 
3.3.2 界面 换 热 系数 的 验证 ”将 式 (20) 和 (21) ( 建立 了 熔 模 铸造 铸件 /型 壳 换 热 一 维 传 热 数 


代入 图 2 所 示 的 铸件 的 三 维 模型 中 , 即 铸 件 表 面 每 学 模型 , 讨论 了 在 导热 系数 和 比热容 不 受 温度 影响 
一 处 的 2 都 根据 表面 单元 及 附近 单元 的 凝固 状态 或 ”的 情况 下 铸 型 内 部 温度 场 与 铸 型 内 壁 热流 密度 之 
表面 单元 的 温度 来 决定 , 其 它 热 物 性 参数 等 均 与 一 间 的 关系 , 从 数学 上 分 析 了 比热容 .导热 系数 热电 
维 反 算 模型 相同 , 可 得 到 铸件 及 铸 型 内 3 根 热电 偶偶 定位 和 热电 偶 测 温 误 差 等 对 界面 热流 和 界面 换 
处 的 温度 变化 . 将 通过 以 上 方法 计算 得 到 的 温度 与 。” 热 系 数 反 算 的 影响 和 它们 之 间 的 不 同 . 

实验 测 得 的 温度 相 比 , 可 得 到 图 14. 从 图 14 可 以 看 (2) 从 数值 模拟 的 角度 评估 了 比热容 .导热 系数 
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和 热电 偶 定 位 误差 对 铸件 、 铸 型 各 处 温度 的 影响 大 
小 , 结合 初步 反 求 结果 对 各 项 误差 进行 了 校正 . 
(3) 获得 了 较为 准确 的 界面 换 热 系数 , 分 析 了 界 
面 换 热 系 数 与 铸件 状态 (凝固 层 厚度 、 温 度 ) 之 间 的 
关系 . 分 析 表 明 , 在 铸件 表面 形成 凝固 层 期 间 , 界面 
换 热 系数 急剧 下 降 约 60%; et 
界面 换 热 系数 下 降 速度 放 组 , 降幅 接近 峰值 的 
20%; 其 后 铸件 已 凝固 , 铸件 温度 下 降 迅 速 , 界面 换 
热 系数 随同 铸件 温度 一 起 下 降 , 但 下 降 速度 比 前 两 
个 阶段 小 得 多 . 
(4) 在 三 维 铸件 / 铸 型 模型 中 对 这 一 关系 进行 了 
验证 . 验证 结果 表明 , 校正 后 的 计算 参数 基本 准确 ， 
反 算得 到 的 界面 换 热 系数 随 铸 件 状态 变化 的 关系 
较为 可 靠 . 
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